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SUrKAitt 


This  report  describes  work  on  two  tasks  related  to  improved 
electrochemical  energy  conversion  systems  for  vehicle  propulsion. 

The  first  is  related  to  the  lithium-air  cell,  specifically  the  re¬ 
duction  of  oxygen  at  a  platinum  electrode  in  1  K  HCIO^  solution  in 
propylene  carbonate.  Voltage  scans  in  the  anodic  direction  at  ?CC  mv/ 
second  gave  an  oxygen  reduction  peak  at  -0,^2  volts  (versus  b.C.E.), 
ar.  anodic  peak  at  -1.3  volts  and  a  shoulder  on  a  cathodic  branch  at 
-1.12  volts  had  not  been  observed  in  the  less  pure  solutions.  There  is 
evidence  of  impurities  even  in  these  highly  purified  solutions. 

The  second  task  reported  on  deals  with  the  mathematical  treatment 
of  battery  data.  The  testing  of  a  battery  to  simulate  an  8-hour 
excursion  at  variovs  load  profiles  is  described.  The  treatment  is 
extended  to  the  simulation  of  vehicle  propulsion  wherein  the  battery  is 
associated  with  a  30  prime  power  source. 


i  - 


Lists  of  Illustrations 


Pege 

Figure  1  Voltage  scens  for  oxygen  reduction  in  propylene 

cerbonete,  seer.  rst*«,  200  mV/sec .  5 

Figure  2.  Load  profile  for  e  battery  opereted  vehicle  during 

6-hour  service .  8 

Figure  3.  Average  load  profile  for  ectual  load  profile  given 

in  figure  2 . 8 

Figure  b.  Actual  load  profile  for  a  50  KW  vehicle  on  8-hour  duty  ..  10 

Figure  5  •  Average  load  profile  for  ectuel  loed  profile  given 

In  figure  U . 10 

Figure  6.  Stepwise  approximation  of  load  profile  of  figure  5  •  •  •  10 

Figure  7.  Probability  graph  of  time  distribution  of  loads  ....  lb 

Figure  8.  Five  loed  profiles  selected  at  random  from  35  sets  of 

random  variables . 15 

Figure  9*  Average  voltage  -time-cycles  surface  .........  17 


ii 


List  of  Tables 


Page 


'fable  1.  Experimental  Boiling  Points  of  Propylene 

Carbonate .  3 

fable  II.  Number  of  profiles  versus  loveet  state  of 

discharge  and  final  state  of  discharge  . .  22 


iii 


I.  IiiTRODUCnOii 


This  is  the.  fifth  semi-annual  report  of  research  on  high-energy 
electrochemical  energy  conversion  systems.  The  overall  program  to  this  point 
has  been  divided  into  seven  tasks: 

1.  betermination  of  open  circuit  potentials  of  a  series  of  couples  in 
various  electrolytes,  and  at  appropriate  temperatures; 

?.  ue termination  of  the  reversibility  of  these  couples; 

3.  Electrochemical  studies  of  high-energy  couples  leading  to  evalua¬ 
tion  of  these  couples  as  materials  for  construction  of  high-energy  electrically 
rechargeable  storage  systems; 

Jetermination  of  kinetic  parameters  and  evaluation  of  the  rate 
limiting  factors  of  selected  reactions  of  electrochemical  couples  at  appropriate 
electrodes; 

5.  ^termination  of  kinetic  parameters  of  selected  reactions  at  catalytic 
electrodes. 

(>.  Investigation  of  ion  transport  processes  in  membranes  and/or  electro¬ 
lytes  at  elevated  temperatures,  and 

7.  iiathematical  analysis  of  performance  characteristics  of  electro¬ 
chemical  energy  conversion  devices. 

The  first  four  semi-annual  progress  reports  described  work  on  Tasks  1, 

2,  3,  ^1  6,  and  7*  buring  this  reporting  period  vork  continued  on  Tasks 
and  7. 
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n.  TASK  FOUR  -  OXICiER  LIECTKOEE  IN  ^RJPTLEKE  CARBONATE 


A.  Introduction 

This  is  e  continuation  of  the  investigation  of  the  oxygen  reaction  at 
a  smooth  pletinum  electrode  in  the  eolution  IK  LiCIO^-  propylene  carbonate. 
The  theoretical  besis  for  this  work  has  been  discussed  (l). 

B .  Sxoe  rlmental 

1.  Purlficetion  of  materials.  Propylene  carbonate  wes  distilled  at 
3  mm  Hg  pressure  in  the  apparatus  previously  described  (2).  This  pressure 
was  selected  to  lower  the  boiling  point  sufficiently  to  avoid  thermal  de¬ 
gradation,  yet  maintain  adequate  boiling  point  difference  between  propylene 
glycol  and  propylene  carbonate  (3).  A  first  cut  amounting  *o  25  percent  of 
the  total  was  discarded.  Boiling  point  data  observed  are  given  in  Table  I. 

Lithium  perchlorete  (u.  r.  Smith  Co.,  trihydrate)  used  in  this 
experiment  was  first  recrystallised  from  weter,  then  heeted  in  vacuo  to 
fO°C  until  the  monohydrete  was  formed.  The  temperature  was  then  reised  to 
220°C  for  several  hours.  The  weight  loss  of  the  resulting  materiel  wes 
sufficient  to  account  for  three  moles  of  weter  being  loet. 

"dltra-high-purity"  nitrogen  and  "aviator's  breathing"  oxygen  (Air 
Products  Comoeny)  were  passed  through  a  drying  tube  containing  phosphorus 
pentoxide  *ust  before  entering  the  cell. 

Solvents  used  in  the  oxygen  determination  experiments  were  reagent 
grade  or  redistilled. 

2.  a  'paratns.  The  apparatus  used  has  been  described  in  previous 
reports  (1,  2). 

3.  Procedure.  The  electrode  was  potentiostated  et  the  initial  volt’ 
and  .as  passed  into  the  solution  through  a  fritted  gless  ges  sparger  for  a 
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TABU)  I 


EXPERIMENTAL  boiling  points  of  propylene  carbonate 


Pressure 

(mm.Hg) 


Heed  Temoereture 
(UC) 


2.0 

3.0 

3.9 

6.5 


76 

8U 

86.5-8? 

93 
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period  of  thirteen  minutes.  The  gae  flow  wee  then  stopped  for  two  minutes  to 
•  How  the  solution  to  come  to  reet  and  the  eean  wae  initiated.  At  least 
three  scans  were  made  in  the  same  manner  before  ’•eproducibility  was  achieved. 
The  first  and  second  scens  of  each  series  always  gave  higher  currents. 

C.  nesults  end  LUecusslon 

beans  were  made  at  100  mv/ second  and  200  mv/second  in  continuance  of 
the  work  reported  in  the  last  report  (3).  Tha  ecens  made  at  100  mv/second 
showed  no  distinct  difference  between  the  oxygen- saturated  solution  and  the 
nitrogen-saturated  solution.  The  scans  made  at  200  mv/see.  are  shown  in 
Figure  1.  An  oxygen  reduction  peak  occurred  at  approximately  -0.67  volt 
versus  the  d.C.E.  When  the  background  was  deducted,  tha  peak  wes  shifted  to 
about  -0.42  volte  and  amounted  to  about  38  microamps/cm^  of  geometric  surfeee. 

It  may  be  observed  in  comparing  this  eean  with  those  previously 
reported  (p)  that  the  currents  obtained  decrease  markedly  as  the  purity  of 
th"  solution  Improves.  The  major  part  of  tha  currents  obtained  were  still 
apparently  the  result  of  impurities  in  the  solution  although  in  this 
experiment  no  peaks  at  -1.8  cathcdic  and  -1,4  anodic  corresponding  to  water 
or  propylene  glycol  respectively  were  present,  A  feature  of  this  scan  not 
observed  in  previous  runs  was  the  shoulder  on  the  cathodic  branch  at  about 
-2.12  volts.  Also  on  this  run  an  anodic  peak  at  -1.3  volts  was  present. 

1‘hese  had  not  been  observed  previously,  but  were  possibly  masked  by  larger 
peaks  from  other  impurities  on  scans  made  in  less  pure  solutions. 

in  work  reported  by  Toni,  et.  al.  (4),  the  oxygen  wave  in  nitro- 
eodimethylamine  eolution  using  phenyltrimethylammonlum  hexafluorophosphate 
as  supporting  electrolyte  occurred  at  -bout  -1.0  volt  compered  with  about 
-0.7  volt  in  this  investigation.  Ka  showri  by  meane  of  the  relationship  of 
peak  current  to  scan  rate  that  the  mechanism  of  the  reaction  in  his  solution 
was  an  ECS  type,  i.e,,  a  chemical  step  between  two  elaetroehemical  steps. 
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100 

50 
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-50 
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Figure  1.  Voltage  scans  for  Oxygen  reduction  in  Propylene 
Carbonate,  scan  rate,  200  aV/sec. 


He  also  suggested,  but  without  supporting  ovidence,  except  that  it  would  be 
expected  to  be  ccnsietent  with  his  resulte,  that  there  existed  a  reaction 
regenerating  0^  from  the  product  of  the  firet  charge-tranefer  step.  The 
amount  of  v-ter  preeent  affected  his  resulte;  increasing  amounts  of  water 
resulted  in  higher  peak  currents  for  the  first  electrochemical  step  and 
less  negative  voltages  for  the  second  cathodic  peak. 

The  major  differences  between  Toni's  solution  and  ours  are  simply 

these : 

a )  different  solvent 

b)  different  electrolyte 

It  would  not  be  expected  that  differencee  in  viecosity  or  solubility 
alone  would  be  responsible  for  our  lower  oxygen  reduction  currents,  for 
Toni  found  a  considerable  portion  of  hie  current  wae  due  to  adsorbed  oxygen 
by  his  experimente  using  different  potentiostatic  holding  times. 

It  muet  be  that  tiiere  exieted  eome  subetance  in  our  eolution  which 
prevented  oxygen  adeorption  onto  the  electrode.  Thie  could  be  one  of  three 
“hings:  propylene  carbonate,  lithium  perchlorate,  or  Impurities.  It  cannot 
be  determined  without  additional  experimental  wo.-k  which  of  these  three 
possibilities  actually  affected  the  results.  This  can  be  checked  by  (l) 
changing  solvents;  (2)  changing  supporting  electrolyte,  and  (3)  improving 
the  purity  of  the  solution  still  further. 

The  current  levele  obtained  by  Toni  (H)  were  on  the  order  of  1  milli- 

2 

amp/cm  .  Since  these  were  peak  voltametric  currents  It  would  be  safe  to 

assume  steady-state  currente  in  the  range  of  two-thirds  as  much,  or  approxi- 
2 

mately  ?00  microamps/cm  . 


m.  TASK  SEVEK  -  KATHEKATICAL  AMUfSIS  OF  ELECTROCHFHICaL 
ENERGY  COMVEKSIOtj  DEVICES 

Battery  Teetlng  to  a  Handcm  Load  Profile 

A.  Cbiectlve 

The  ideal  evaluetion  of  a  battery  for  e  specific  eppli cation  is  to 
test  the  battery  under  conditions  identical  to  those  under  which  it  will  be 
used.  However  the  loed  profile  of  ectual  use  conditions  raey  be  difficult  to 
determine.  For  example,  for  e  battery  opereted  vehicle  driving  over  en 
"average"  terrain,  e  power  vs.  time  curve  would  be  neerly  random  in  power 
requirements  due  to  accelerations,  stops,  etc.,  and  would  be  difficult  to 
c.  scribe.  This  real  time  load  profile  could  be  in  the  form  of  Fig.  2  over 
an  8-hour  driving  period. 

Over  a  period  of  time,  say  8  hours,  en  estimate  cen  be  made  of  the 
total  time  spent  et  e  given  loed  end  an  everage  loed  profile  determined. 

For  the  vehicle  profile  of  Fig.  2  ,  the  average  load  profile  fc  8  hours  would 
be  of  a  form  similar  to  Fig.  3  where  certain  loads  ere  more  frequent  than 
others.  A  loed  profile  for  battery  testing  cen  be  created  by  selecting  loads 
from  the  average  profile  end  applying  them  to  the  bettery  for  the  indicated 
length  of  time ,  The  results  of  such  e  test  cen  be  considered  es  an  indica¬ 
tion  of  how  the  bettery  would  heve  responded  in  the  reel  time  profile,  or  if 
severe  loads  ere  selected,  conclude  that  any  bettery  performing  satisfactorily 
in  this  test  would  elso  perform  satisfactorily  in  reel  time  profile.  Which¬ 
ever  interpretation  is  used  for  the  test  resulte,  the  results  will  be  effected 
by  the  length  of  time  each  loed  is  applied  end  by  the  sequence  of  loads, 
e.g,,  the  higher  loads  in  the  beginning  or  end  of  the  cycle. 

The  test  results  of  e  bettery  will  be  biesed  by  the  particular  sequence 
of  loads  selected,  perticulerly  if  meny  cycles  of  the  seme  sequence  are 
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Figure  2,  Load  profile  for  a  battery  operated 
vehicle  during  8-hour  service. 


KW 


Figure  3*  Average  load  profile  for  actual 
load  profile  given  in  figure  2. 


applied.  But  this  load  sequence  is  abstracted  from  en  average  load  profile 
and  there  are  no  criteria  available  for  selecting  one  sequence  over  another. 

For  instance,  there  is  no  justification  for  selecting  a  test  profile  composed 
of  a  high  power  step  followed  by  a  low  power  step  followed  by  a  medium  power 
step.  Likewise  the  length  of  time  which  e  given  load  is  applied  is  an 
estimate  of  an  average  and  need  not  be  e  fixed  number.  Since  there  is  no 
reason  for  selecting  any  given  sequence,  the  most  meaningful  and  useful  date 
can  be  obtained  by  testing  the  battery  under  a  variety  of  load  sequences 
and  durations  of  load  application.  It  is  possible  to  simulate  this  variety 
by  Lonte  Carlo  methods. 

B.  Discussion 

A  load  profile  and  S-hour  discharge  sequence  postulated  for  a  50  W 
vehicle  which  corresponds  to  the  equivalent  of  four  full  power  houra  isi  (5) 

a)  50  KW  for  1  hr 

b)  25  KW  for  1  hr 

c)  75  KW  for  1/3  hr 

d)  1*5  KW  for  2/3  hr 

e)  10  KW  for  2  hrs 

f)  Open  Circuit  for  2  hrs 

g)  50  KW  for  1  hr 

Figure  4  is  a  graphical  representation  of  this  profile.  However,  repeated 
cycling  of  a  battery  to  the  same  load  profile  may  produce  results  unique  to 
that  particular  test  profile,  and  not  a  good  indication  of  how  the  battery 
would  behave  in  actual  use.  A  variety  of  discharge  sequences  is  needed  to 
eliminate  test  artifacts  peculiar  to  e  given  sequence. 

The  above  load  profile  (Fig.  4)  can  be  mede  more  general  by  considering 
it  to  have  been  formed  from  the  8-hour  average  load  profile  of  Fig.  5  which 
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Figure  4.  ftctuai 
load  profile  for 
a  50  vehicle  on 
fi-hour  duty 


Figure  5*  Average 
load  profile  for 
actual  load  profile 
given  in  figure  h 


Mgure  6.  itepvdfe 
approximation  of 
load  profile  of 
figure  5 


ie  the  data  of  reference  5  with  the  set  sequsnos  eliminated  a«i  expressed 
as  a  continuous  curve. 

Fig.  5  can  be  step-wise  apprordnabed  into  Fig.  6. 

,is 

fhis/now  in  the  form  of  an  8-hour  average  load  profile  similar  to 
Fig,  3  and  new  discharge  eequencse  can  be  selected.  To  generate  a  naw  se¬ 
quence,  the  first  load  to  be  applied,  the  escond  load,  etc.  must  be  selected 
by  some  method.  Since  there  is  no  juetification  for  selecting  sny  load 
preferentially,  loads  must  be  selected  st  random.  Then,  if  the  numbers 
/0,  1,  2,  ---,  2 7  are  assignad  to  correspond  to  the  load*  ,/OKW,  1C  KW,  20  KW,--- 
70  KW/,  a  random  selection  of  a  number  Jp,  1,  2  will  determine  the 

first  load  of  the  discharge  sequence.  Similarly,  a  second  random  number  is 
salacted  to  determine  the  second  load  of  the  eequance,  etc.  If  each  load 
Is  to  occur  only  once  in  an  8-hour  discharge  cycle,  the  number  corresponding 
to  that  load  must  be  excluded  from  succeeding  seleotions  once  it  hae  been 
chosan.  A  profile  is  completed  when  all  eight  loads  have  been  applied.  Thus, 
numerous  discharge  sequence s  csn  be  generated  which,  if  applied  succeseivaly 
to  a  battery,  eliminate  sny  bise  to  the  results  due  to  any  one  discharge 
sequence. 

Further  generality  can  be  added  to  the  test  profile  if  the  length  of  time 
which  a  given  load  is  applied  is  not  held  ae  a  fixed  number.  The  test  profile 
is  obtained  from  an  average  load  profile  which  ie  an  estimate  of  the  langth 
of  tima  which  a  given  load  is  applied  ovar  an  8-hour  period.  This  length  of 
time  need  not  be  the  same  for  avery  8-hour  period  but  can  be  considered  to 
be  an  average  value  of  many  8-hour  periods.  Than,  if  many  8-hour  periods  ara 
taken,  the  langth  of  time  a  given  load  is  applied  will  be  a  variable  which 
has  an  average  and  is  distributed  about  this  avsraga  by  some  probability 
distribution.  This  will  allow  a  given  load  to  be  applied  for  numerous  time 
durations  which  will  tend  toward  an  average  value. 
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For  convenience,  essume  this  probability  distribution  to  be  Normal 
(Geusslen).  An  average  end  a  standerd  deviation  (o)  are  neceesary  to 
describe  the  tine  allotted  for  each  of  the  eight  loade  mentioned  above.  The 
average  tine  will  be  taken  es  that  time  for  each  load  in  Fig.  6.  ^ince  e 
standerd  devietion  is  not  apparent  from  the  average  load  profile,  it  mey  be 
arbitrarily  assigned. 

The  selection  of  o  muet  be  such  that  reasonable  tine  durations  ere 
obtained.  For  Instance,  if  the  average  time  for  a  given  load  Is  2  hours  it 
would  not  be  expected  that  in  eny  one  8-hour  profile  this  tine  would  be 
drastically  different  from  2  hours,  but  would  usually  be  between  1.5  hours 
to  2.5  hours.  If  o  is  chosen  to  be  one-eighth  the  average  tine  for  each 
load,  reasonable  time  distributions  ere  obtained. 

Then,  by  the  properties  of  the  Normal  distribution,  in  96 $  of  the 
8-hour  average  load  profiles,  the  length  of  time  which  e  given  load  is 
applied  will  be  within  /Avg.  1  2o_/.  For  the  eight  loads  discussed  above, 
the  tine  values  will  be: 


Loed  (KW) 
0 


/Avg.  -  ZqJ  (h  rs ) 


2  -  .5 


10 

1.3?  !  03 

20 

l.o  1  .25 

30 

•33  !  .08 

40 

.6?  1  .1? 

50 

2  !  .5 

60 

•33  -  .08 

70 

•33  -  .08 
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Fig.  7  Is  a  probability  graph  of  the  tins  distribution  for  these 
loads.  It  is  read  as:  for  any  load,  what  is  the  probability  that  the  load 
will  be  applied  for  a  length  of  tine  leea  than  or  equal  to  aone  selected  tine. 
Then  a  length  of  tine  to  apply  a  load  will  be  determined  by  randomly  selecting 
a  two  digit  number  £, 00,  .01,  .02,  — ,  1.0 ]  to  determine  a  probability  value 
and  read  from  Fig.  7  the  length  of  time  corresponding  to  that  particular 
load.  For  instance,  if  the  number  selected  to  correspond  to  the  load  is  5 
(50  KW)  and  the  probability  value  selected  ie  .25,  a  50  KW  load  ia  applied 
for  1.84  hr  as  the  first  step  of  a  discharge  sequence.  A  new  load  and 
probability  value  is  selected  for  the  second  step  of  the  sequence,  etc. 

Thirty-five  sets  of  random  variables  were  chosen  to  determine  five 
different  test  load  profiles  and  are  shown  in  Fig.  8.  There  is  no  noticeable 
pattern  of  occurrence  of  any  particular  load  among  the  cycles. 

The  length  of  any  cycle  should  be  approximately  eight  hours.  The  time 
limits  of  the  generated  cycles  are  7*12  hra.  and  8.36  hra.,  the  average  of 
the  five  is  ?»97  hrs.  The  average  length  of  time  each  load  is  applied  is: 

Load  (ICiQ  Avg,  Time  (rirs) 


0 

1.97 

10 

1.40 

20 

1.01 

30 

0.32 

40 

0.62 

50 

2.09 

60 

0.31 

70 

0.32 

wiich  are  approximately  the  values  of  Fig,  6, 
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Each  of  these  test  load  profilee  removes  approximately  the  seme 
emount  of  power  (everage  of  the  five  being  205.1  watt-hrs)  but  in  e 
different  sequence  of  dischargee. 

It  ie  neeeesary  to  collect  the  teat  result*  into  a  for*  which  indicate e 
battery  reeponee  to  the  average  load  profile  (Fig.  6)  but  without  a  parti¬ 
cular  discharge  eequsnce.  l«et  us  assume  it  ie  neeeesary  to  evaluate  battery 
response  to  a  y>  KW  load.  Thle  responee  will  dspend  on  the  loads  prior  to 
the  50  KW  load  and  on  the  number  of  cycle*  on  the  battery.  Since  eny  number 
and  velue  of  loade  could  precede  the  50  KW  load,  all  possible  prior  loede 
muet  be  considered.  Thle  is  the  extremely  large  number  of  ell  combinations 
and  permutations  of  the  other  seven  load  values.  The  battery  responee  to 
a  5°  KW  load  will  then  bs  the  average  of  the  response  to  all  poesible  prior 
loads. 

If  the  five  test  load  profilae  of  Fig.  8  are  used,  an  approximation 
to  all  poeelble  loads  prior  to  50  KW  ie  made.  If  the  voltage- time  curves 
of  the  five  load  profiles  are  obtained,  the  bet.t  estimate  of  the  voltage 
during  a  50  KW  losa  will  be  the  average  voltage  for  the  50  KW  load  of  ths 
five  curves.  This  seme  averaging  technique  is  used  on  all  eight  load  values 
to  give  a  voltage- time  curve  after  five  cycles  which  ie  independent  of 
discharge  sequence,  and  can  be  arranged  into  a  voltage  reeponee  curve  corres¬ 
ponding  to  the  power  requirements  of  Fig.  6. 

After  numerous  cycles  sre  applied  to  the  battery,  the  data  is  collected 
in  five  cycle  sets,  and  averages  of  the  responses  to  each  load  are  taken. 

This  represents  the  battery  response  between,  for  example,  cycles  51  to  55. 

The  same  procedure  is  continued  to  compile  a  voltage -time- cycles  contour 
similar  to  Fig.  9* 
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The  experimental  reeults  of  battery  testing  are  now  in  e  form  in 
which  the  results  are  independent  of  the  discharge  sequence  and  present  e 
more  reliable  picture  of  the  effect  of  bettery  cycling. 

Because  of  the  symmetry  of  voltage-time  end  curren'Utime  curves 
during  constant  power  discharges,  e  current- time-cycles  contour  can  be  eesily 
constructed  from  the  voltage  contour. 

Vehicle  Simulation 

A.  Discussion 

Thus  fer  it  has  been  implied  that  the  power  for  vehicle  propulsion 
will  be  supplied  solely  by  batteries  which  would  be  recherged  at  the  end 
■  f  the  8-hour  driving  period.  In  a  hybrid  system,  e  prime  power  source  is 
provided  in  addition  to  batteries,  such  thet  battery  power  is  required  only 
for  the  higher  power  demands  of  the  vehicle.  At  low  vehicle  power  demands, 
power  from  the  prime  source  not  required  for  vehicle  propulsion  is  eveilable 
for  bettery  recharging.  A  rendom  sequence  of  the  load  profile  of  Fig.  8 
epplied  to  a  hybrid  system  would  then  ceuse  e  rendom  mixture  of  charges 
end  discherges. 

It  is  desired  to  select  a  prime  power  source  and  e  secondary  bettery 
S’.ch  thet  sufficient  power  is  available  to  meet  the  vehicle  load  profile  for 
all,  or  nearly  ell  loeds.  This  selection  ie  subject  to  minimizing  the  power 
requirements  of  the  bettery  for  a  given  prime  power  source  and  requiring 
that  the  better;-  be  recharged  in  e  reasoneble  length  of  time  et  the  end  of  a 
driving  period. 
b.  .-tesults 

A  50  KW  vehicle  is  essumed  and  the  load  profile  of  the  vehicle  is 
that  of  rig.  6.  It  is  necessary  to  select  a  prime  power  source  and  secondary 
battery  combination  which  meet  the  criterie  described  ebove.  As  e  first 


estimation,  a  prime  power  source  of  30  KW  is  assumed.  It  is  then  possible 
to  calculate  what  battery  size  would  be  optimal  with  this  size  prime  power 
source,  .tith  the  30  hi.  prime  source,  discharge  loads  of  thp  battery  are 
then  iiO  iC.V,  30  K.»,  20  KW,  10  KW;  battery  charging  levels  are  30  KW,  20  ICW, 

10  KW.  At  the  30  KW  vehicle  load,  the  battery  is  neither  charging  nor  dis¬ 
charging.  a  ..ante  Carlo  simulation  of  vehicle  operation  was  made  by  forming 
ZC  random  sequences  of  the  eight  level  load  profile.  Calculations  were 
performed  on  a  PDP-8  computer  to  determine  the  KWH  of  charge  or  discharge 
applied  to  or  removed  from  the  battery  by  each  step  in  the  sequence,  and 
to  determine  by  a  cumulative  sum  the  state  of  discharge  of  the  battery  during 
the  sequence. 

Topical  of  these  calculations  is  as  follows: 


!iicle 

}  c  it*.*  J 

battery* 
Load  (KW) 

Time  (Hrs) 

KWH 

Ea  tte  ry  * 
Citate  (KWH) 

10 

+20 

1.33 

00 

00 

00 

+30 

2.00 

00 

00 

5C 

-20 

2.00 

O 

1 

-L0.0 

70 

-90 

0.33 

-13.3 

-33.3 

20 

+10 

1.00 

+10 

-*3-3 

60 

-30 

0.33 

-10 

-33-3 

30 

00 

0.33 

00 

-53.3 

if  0 

-10 

0.6? 

-6.7 

-60.0 

*  +  =  charge,  -  =  discharge 

In  the  particular  profile,  the  first  vehicle  load  is  10  KW,  but  since 
30  K..  is  available  from  the  urine  source,  2C  K«  of  charge  is  applied  to  the 
battery  for  1.33  hours.  However,  as  the  battery  is  assumed  fully  charged 
at  the  beginning  of  the  day,  no  charge  is  given  to  the  battery,  so  the 
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battery  state  is  zero  (fully  charged).  In  the  third  step  of  the  profile 
50  Kn  is  demanded  by  the  vehicle  eo  20  K M  must  be  removed  from  the  battery' 
for  2  hours  which  removes  40  KWH  from  the  battery  in  this  step,  leaving  the 
battery  state  40  KWH  discharged.  Similarly,  the  rest  of  the  table  is 
interpreted. 

Tne  data  from  each  of  the  56  profiles  were  examined  for  the  lowest 
state  of  discharge  in  each  sequence  and  the  stete  of  discharge  at  the 
completion  of  the  6-hour  period.  The  number  of  profiles  having  a  given 
lowest  stete  of  battery  discharge  and  the  final  state  of  discharge  is 
tabulated  (expressed  as  probability)  in  Table  II. 

Irom  the  teble  it  is  apparent  that  the  battery  is  moet  likely  to  end 
the  2-hour  period  in  a  fully  charged  state.  The  next  most  likely  final 
state  of  the  battery  is  a  40-50  KWH  discharge  6tete.  The  40-50  KWH  dis¬ 
charge  state  is  also  the  most  likely  deepest  stete  of  discharge;  however,  a 
60  iCi.H  or  greater  stete  of  discharge  occurs  a  significant  number  of  times. 

The  high  stete  of  discharge  of  the  battery  is  less  importent  if  the 
battery  does  not  reach  that  stete  until  the  very  end  of  the  profile,  as 
the  battery  will  be  recharged  regardless.  However,  from  the  chart  it  is 
apparent  that  the  likelihood  of  not  reaching  the  lowest  state  until  the  end 
of  the  profile  is  quite  small. 

In  73  of  the  86  profiles,  a  "waste  of  charge"  was  involved,  i.e.,  a 
charge  could  have  been  epp’ied  to  the  battery  et  some  step  in  the  profile 
if  the  battery  were  not  already  fully  charged.  Sixty-six  of  these  73  profiles 
involved  a  "waste  of  charge"  at  other  than  the  laet  step  of  the  profile.  The 
distribution  of  the  amount  of  charge  wasted  during  these  66  profiles  is: 
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Amount  of  charge 

"wasted"  (KWH)  A.  ?0  70-60  60-50  50-1*0  40-30  >0-20  20-10  10-0 

Fraction  of  profiles 
wasting  charge  of 

amount:  .06  .05  .23  .21  .17  .11  .11*  .03 

If  this  "wasted"  charge  could  have  been  utllliad,  e.g.,  having  an 
auxiliary  battery  which  could  always  be  charged,  the  lowest  state  of  discharge 
of  the  battery  would  have  been: 


State  of  discharge 

(KWH)  ±60  60-50  50-1*0  1*0-30  30-20  20-10  10-0  Charged 


Fraction  of  profiles 
whose  lowest  state 

of  discharge  is:  .14  .23  .29  .09  .05  .15  *05  0 
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TABLE  II 


FHACTION  OF  PROFILES  VERSUS  LOWEST  STATE 
OF  DISCHARGE  ADD  FINAL  STATE  OF  DISCHARGE 


State  of 
discharee  (KWH } 

Fraction  of 
profilss  with 
final  state 
of  discharee 

Fraction  of 
profilss  with 
lowest  state 
of  discharee 

Fraction  of 
profiles  with 
lowest  state 
being  final 
state 

ioO 

.05 

.17 

.05 

60-50 

.04 

.33 

.09 

50-40 

.19 

.46 

.16 

40-30 

.0? 

0 

0 

30-20 

•  0? 

.01 

0 

20-10 

.16 

.01 

0 

10-0 

.09 

0 

0 

Charged 

.26 

0 

0 
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Apparently  some  improvement  is  made  in  reducing  the  lowest  state 
of  discharge,  e.g.,  the  probability  of  the  lowest  discharge  being  greater 
than  50KW  is  decreased  from  .50  to  .37;  however,  a  significant  probability 
exists  in  having  a  discharge  otate  groater  than  60KW. 

C.  Conclusions 

Simulation  Involving  86  load  profiles  of  the  81  possible  load  profiles 
does  not  ellovi  other  than  the  roughest  conclusions  to  be  drawn.  From  this 
simulation  which  assumed  a  30KW  prime  power  source,  it  appears  that  a 
battery  capable  of  approximately  60  KWH  is  necessary  to  meet  most  vehicle 
load  profile  requirements  during  the  8-hour  drive.  Further,  at  the  end  of 
the  8-hour  period,  the  battery  Is  likely  to  be  in  e  charged  state  or  in  a 
discharge  state  of  less  than  approximately  hO  KWH, 

While  increasing  the  power  of  the  prime  source  will  aid  in  reducing 
the  battery  power  requirements,  it  will  also  increase  the  likelihood  of  the 
battery  ending  the  profile  in  e  fully  charged  state,  Implying  Inefficient 
use  cf  the  battery.  It  appears  then  that  an  optimum  selection  of  prime  power 
source  and  battery  power  must  exist  but  further  constraints  such  as  cost, 
weight,  etc.  must  be  included.  Further  refinements  in  the  simulation 
calculations  end  inclusion  of  additional  constraints  is  not  likely  warranted 
until  a  real-time  load  profile  becomes  available. 

In  a  similar  manner  to  that  eerlier  descrioed,  five  load  profiles  in 
which  a  hybrid  configuration  is  assumed  can  be  selected  for  battery  testing, 
rron  the  test  results  voltage-time-cycles  contours  will  be  obtained  but 
in  addition,  overcharge  characteristics  of  the  battery  can  be  evaluated. 
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